Thermally activated switching in the presence of non-Gaussian noise 
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We study the effect of a non-Gaussian noise on interstate switching activated primarily by Gaus- 
sian noise. Even weak non-Gaussian noise can strongly change the switching rate. The effect is 
determined by all moments of the noise distribution. The explicit analytical results are compared 
with the results of simulations for an overdamped system driven by white Gaussian noise and a 
Poisson noise. Switching induced by a purely Poisson noise is also discussed. 
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Much progress has been made recently in the studies 
of switching between coexisting stable states, primarily 
because switching can be now investigated for a large 
variety of well-controlled micro- and mesoscopic systems 
ranging from trapped electrons and atoms to Joscphson 
junctions and to nano- and micro-mechanical oscillators 
[U, i, i, i, i, i, 0, S, H, ■ Fluctuations in these systems 
are usually due to thermal or externally applied Gaus- 
sian noise. However, as the systems become smaller, an 
increasingly important role may be played also by non- 
Gaussian noise. It may come, for example, from one or a 
few two-state fluctuators hopping at random between the 
states, in which case the noise may be often described as 
a telegraph noise. 

The switching probability is sensitive to a non- 
Gaussian noise. This sensitivity attracted much at- 
tention after it was proposed to use switching in 
Josephson junctions to measure the full counting statis- 
tics in electronic circuits [H], [Tsf. Several theoretical 
[13, [H, [H, [l3 ^^'^ experimental [l8l. [l^ papers on mea- 
suring the 3rd moment of the current distribution from 
the switching rates were published recently, and different 
theoretical approaches were compared in Refs. [2Q, 21. 

In this paper we study switching induced by Gaus- 
sian noise in the presence of an additional non-Gaussian 
noise. Even where the latter has a smaller intensity than 
the Gaussian noise, its effect on the switching rate may 
be exponentially strong. We show that it can be de- 
scribed in a simple form in terms of the noise character- 
istic functional, thus accounting for all moments of the 
noise distribution. The analytical results are compared 
with simulations for an overdamped system driven by 
white Gaussian noise and a Poisson noise. We also con- 
sider switching induced by a Poisson noise alone; here, 
the result for the rate may be qualitatively different from 
that in the weakly non-Gaussian noise approximation. 

The potentially strong effect of an extra modulation, 
whether random or regular, on the rate of Gaussian 



noise induced switching can be understood from the well- 
established picture of the switching dynamics. Switching 
events result from large rare noise outbursts. For Gaus- 
sian noise, the switching rate is W oc exp{—R/D), where 
R is the activation energy and D is the noise intensity 
[2^ . Even though switching happens at random, the 
system trajectories followed in switching form a narrow 
tube in the space of dynamical variables q = {qi,q2, ■ ■ ■) 
centered at the most probable (optimal) switching path 

qopt(0- 

One can think of the effect of an additional modulation 
in terms of a generalized work done by the modulation 
on the system moving along qopt {t) [l^, US] • This work 
changes the activation barrier. The change SR is propor- 
tional to the modulation amplitude. Therefore the over- 
all change of the switching rate oc cxp{—SR/D) depends 
on the modulation amplitude exponentially. The switch- 
ing ra te g ives the probability current from the occupied 
state |25| . making it an observable quantity. Thus, if 
the modulation is random, one has to simply average the 
factor exp{—SR/D) over realizations of the modulation. 

Since SR is linear in the characteristic amplitude of 
random modulation, the ratio SR/D does not have to 
be small even where the modulation intensity, which is 
quadratic in the amplitude, is smaller than D. However, 
the distribution of non-Gaussian modulation may decay 
slower than Gaussian on the tail. To determine whether 
the effect of a non-Gaussian noise on switching may be 
regarded as a perturbation one has to compare the prob- 
abilities of appropriate large fluctuations induced by the 
Gaussian and non-Gaussian noises, taking into account 
all moments of the distribution. 

We study switching for a system described by the 
Langevin equation 



q = K(q)+f(i)+|(t). 



(1) 
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We assume that, in the absence of noise, the system has 
a stable stationary state q^A and a saddle point q^ on the 
boundary of the basin of attraction to q^, with K(qyi) = 
K(q5) = 0. Switching from the stable state is due to the 
forces f(t) and ^{t), which are the Gaussian and non- 
Gaussian noises, respectively. We separate them, since 
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physically they often come from different sources. It is 
convenient to characterize f (t) by its probability density 
functional 7'f[f(0] exp (-7^f /£>), 



7^f[f(t)] = \ 



dtdt'f{t)T{t-t')i{t'), 



(2) 



where jF{t — t')/2D is the inverse of the pair correlator 
of f{t). The characteristic noise intensity D is small, so 
that the switching rate W ^ ,t~^ , where is the 
relaxation time of the system and tc is the noise correla- 
tion time. The non-Gaussian noise is more conveniently 
described for our purpose by the characteristic functional 



T'lM = (exp 



dt\i{t)$,{t) 



(3) 



means averaging over ^(t). 



where (. 

We first consider the case where the intensity of the 
non-Gaussian noise £ is smaller than D. We will disre- 
gard corrections proportional to this intensity, but the 
ratio of the characteristic amplitude 50 of ^ to I? will not 
be assumed small. The switching rate can be written as 

W = C{eM-m]/D])i. (4) 
n[i] = min |7^f + i j dtk{t) [q - K - f(t) - , 

where C is the prefactor that weakly depends on the 
noise intensity. The minimum is taken over trajec- 
tories f (t), q(t), k(t) that satisfy boundary conditions 
{{t),k{t) ^ ior t ^ ±00, qt-,_oo <U,<it~,oa ^ <ls- 
This formulation was proposed in the weak-noise limit 
j26l [27| for a time-periodic ^(i), in which case q_A,s are 
also periodic and there is no averaging over ^. The varia- 
tional problem ^ describes coupled optimal trajectories 
fopt(i), qopt(i), kopt(i)i with fopt(i) being the most proba- 
ble noise realization that brings the system to the saddle 
on the basin boundary of the initially occupied state. 

It is known from variational calculus that, to first order 
in ^, the effect of ^{t) on can be calculated along the 
optimal trajectory unperturbed by S,{t). Such a trajec- 
tory is an instanton. Its typical duration is ~ max(<r, tc). 
It is translation-invariant with respect to time and can be 
centered at any time tp. If ^(t) is periodic, it lifts time- 
translation symmetry and fixes tg (modulo the period) 
so as to minimize 7?,[^]. 

If ^{t) is a stationary noise, the switching rate W is in- 
dependent of time. In this case one can think not of the 
adjustment of the instanton center to to but, equiv- 
alently, of the adjustment of ^(t) to to so as to maximize 
the overall probability of switching. This adjustment pro- 
vides the major contribution to the value of W when the 
averaging over realizations of ^{t) is performed in Eq. (P]) 
using a solution with a given to- 

From Eqs. ([4]) one obtains a simple expression for 
the switching rate 



where is the switching rate in the absence of non- 

Gaussian noise. The factor Agv, describes the effect of 
non-Gaussian noise. It is expressed in a closed form in 
terms of the noise characteristic functional calculated for 
function x{t) = ~*kQpj(i), where k.l^J^{t) is the solution of 
the variational problem ^ for ^ = 0. The real function 
x{t) is the logarithmic susceptibility which describes the 
linear response of log W to a perturbation [H, [U, HI] . 
The structure of Eq. ^ resembles that of the expression 
for a large fluctuation probability in a birth-death system 
with non-Gaussian modulation of reaction rates [29|. 

From Eq. the effect of a non-Gaussian noise on 
the switching rate is determined by the ratio of the noise 
amplitude to the Gaussian noise intensity D. Equa- 
tion ([5]) applies to both underdamped and overdamped 
systems. Examples of calculating k^pj can be found in 
Refs. [2^ 113, m and papers cited therein. 

As an illustration we will consider the case of a one- 
component 5-correlated Poisson noise ^ {t) with pulse area 
g and mean frequency v. Using the explicit form of the 
noise characteristic functional l30l. we obtain 



A. 



cxp 



/ dt [I -c^j> {-x{t)g/D)] 



(6) 



A^^^r^iix/D], 



(5) 



where x(^) is the corresponding component of the loga- 
rithmic susceptibihty. If g/D <C 1, logAgw is a series in 
g/D. The coefficients in this series describe the effects 
of the moments of the Poisson noise on the switching 
rate. In the opposite case, 5/D ^ 1 (but the Poisson 

noise intensity vg^ <C D), if xit) becomes negative, then 

1/2 

\ogAs„ « z/ [27r£)/gx(im)] exp[-x{tm)g/D] where t^ 
is the instant where — x(*) is maximal. If x{t) ^ for 
all t and g/D 3> 1, the major contribution to Agw comes 
from the region of small x(^)- If x(0 is small only for 
\t\ — > 00, where it decays exponentially with \t\, then 
log oc i/tr \og{g/D), to leading order in g/D. 

The Poisson noise distribution does not fall off as 
steeply as Gaussian. This imposes a limitation on the 
range of g/D where Poisson noise may be treated as a 
perturbation and the above theory applies. To see the 
far-tail effect we consider switching due to a purely Pois- 
son noise, where f = in equation of motion ([1]). We 
will use the method of optimal fluctuation, as for some 
other types of non-Gaussian noise (3l| . 

The switching rate is determined, to logarithmic ac- 
curacy, by the integral over trajectories k(t),q(t) of the 
functional (exp{-i/dtk(t)[q-K-|(t)]})^ [3|. We 
assume for brevity that different components of the Pois- 
son noise are independent short pulses with areas g = 
(gi, 52, • ■ •) and with average frequency v. It is convenient 
to consider a zero-mean noise, (,{t) ^{t) — lyg, K — > 
K + i^g; we assume that qA,q5 are also appropriately 
shifted. Of interest for switching are trajectories that ap- 
proach the saddle point [1^ . In the spirit of the method 
of optimal fluctuation, for small |g| and for i' t~^ the 
integral over trajectories q(t),k(t) can be calculated by 
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steepest descent. This gives 

W = C"exp[-i?p], Rp = min J dt (ikq - H) ,(7) 

H = H{q, ik) = -u{l + igk - e'sk) + ikK(q). 

The variational problem ([7]) determines the optimal 
switching trajectory qopt(Oi ^pt{t)- It starts at t —f —oo 
at q ^ q^ijk ^ and goes to q ^ q5,k ^ for 
t — *■ c». On this trajectory H = 0. As in systems 
driven by white Gaussian noise [l^l, the optimal tra- 
jectory can be described as a Hamiltonian trajectory of 
an auxiliary system with coordinate q, momentum ik, 
and Hamiltonian H. A similar approach was proposed 
m Ref. [H. However, the characteristic functional was 
not specified and the explicit analysis took into account 
only the 3rd moment of the noise distribution, which was 
considered a perturbation; therefore the results do not 
describe switching due to large Poisson fluctuations. 

The switching exponent in Eq. ([7]) is i?p ^ 1 for small 
|g|. However, in contrast to the case of Gaussian noise, 
Rp is not proportional to the reciprocal noise intensity 
i/g^. Nor does it scale like reciprocal noise amplitude 
|g|~^, although i?p |g| appears to slowly vary with |g|. 

An explicit dependence of the switching rate on the in- 
tensity g of Poisson-distributed pulses can be found for a 
one- variable overdamped system with equation of motion 



'U\q) + .f(t)+at). 



(8) 



Here, U{q) is the effective potential. The stationary 
states qA and qs correspond to the minimum and the 
barrier top of U{q). 

If f{t) is white Gaussian noise, {f{t)f{t')) ^ 2DS{t - 
t'), and the Poisson noise is weak, the Poisson-noise in- 
duced factor in the switching rate is described by 

Eq. © with x(t) = -/ipi(i)/2 = aiid with 

In the opposite case where switching is due to purely 
Poisson noise, i.e., / = in Eq. (O, from Eq. ([7]) 



.(0) 
Vopt 



Rp 



dqn{q), 



\og{l + [KU'{q) I gi^]]. (9) 



Here, qA and qs are the shifted extrema of the potential 
given by equation U'{q) = gv. From Eq. ([5]), i?p ~ 
rp \og{rp/vt^), with rp = {qs - qA)/g- 

A qualitative feature of unipolar (pulses of one sign) 
Poisson noise is that, for an overdamped system, it causes 
switching only provided the noise pulses push the system 
from the stable state towards the saddle. In this case 
rp > 0. There is no switching for pulses of the oppo- 
site sign. The "one-sidedness" of fluctuations in over- 
damped systems has other manifestations, which includes 
the work fluctuation distribution [3^ . On the other hand, 
we expect that an undcrdamped system should be able to 
switch for Poisson pulses of any sign, in which case there 
should be a critical value of damping for which switching 



from the state is possible for a given sign of g. Equations 
©, ® apply iirp,rp/vt, > 1. 

As a cause of switching, a Poisson noise is effec- 
tively weaker than a Gaussian noise if the switching 
exponent Rp is larger than the switching exponent for 
the Gaussian noise. For white Gaussian noise of in- 
tensity D in Eq. ([5]), the switching exponent is ISU/D 
with MJ = U{qs) - U{qA) m. The condition Rp > 
AU/D effectively limits the range of applicability of 
Eq. We saw that logAgw becomes large provided 

—x{tm)g/ D ^ 1. An order of magnitude estimate shows 
that ~x{tm)g/D - {AU/D)Rp^\og{rp/iytr), and there- 
fore from Eq. ^ the large log^sw asymptotics applies 
only provided \og{rp/i>tr) 3> 1. This condition is com- 
patible with —x{tm)g/D ^ 1 only for very small D. 

We now apply the above results to an overdamped sys- 
tem dH]) with a double-well potential 



(10) 



which has been extensively studied in the context of 
white-noise driven systems. In the absence of Pois- 
son noise, the escape rate in this case is W^^"^ ~ 
(v^/tt) exp(— 1/41?), and the logarithmic susceptibil- 
ity for escape from the negative-q well {qA = —1) is 
x{t) ~ — exp(i/2)(2coshi)~'^/^ (x(0 opposite sign 
for switching from qA = 1). 
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FIG. 1: Poisson noise induced change of the switching expo- 
nent for an overdamped Brownian particle in potential (|10|l . 
The dimensionless Poisson noise frequency is v — 0.5. The 
data of numerical simulations are shown with circles and 
squares for the cases where in escape the particle moves along 
and opposite to Poisson pulses, respectively. The solid curves 
show the weak-Poisson noise theory ([6]) for these cases, and 
the dashed curves show the approximation where only the 
three moments of Poisson noise are taken into account. 

In Fig. [1] we present results of Monte Carlo simula- 
tions of switching of an overdamped Brownian particle 
described by Eqs. ([8]), (fTO|) . They are compared with 
the weak Poisson noise prediction, Eq. ([5]), and with the 
approximation where only terms up to g^ are kept in 
Eq. ([6]). The Poisson noise intensity I'g^ < D in the 
whole range of studied g/D. For g/D < 1 the effect of 
Poisson noise is small and comes, primarily, to the change 
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of the activation barrier, AC/ — s- AC/ ± vg, and the effec- 
tive noise intensity, D — > D+vg^ . For larger g/D > 3 the 
switching exponent changes significantly, as expected. 

For switching in the direction opposite to Poisson noise 
pulses, where \ogW/W^^^ < 0, the numerics agrees well 
with Eq. This is to be expected, since the far 

tail of Poisson noise distribution is immaterial here; for 
g/D > 5 the results differ noticeably from the three- 
moments approximation. For switching along the pulses, 
because of the far-tail effect, with increasing g Poisson 
noise quickly becomes as important as white noise for 
chosen D. Therefore the perturbation theory fails and 
the dependence of the switching exponent on g is much 
weaker than the exponential dependence expected from 
Eq. (El). 

20 1 ' ' ' ' 1 
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FIG. 2: Switching exponent Rp = —InW for an overdamped 
particle in a potential (|10|) driven by a Poisson noise with 
dimensionless mean frequency v = 0.5 and pulse intensity g. 
Squares show the results of Monte-Carlo simulations, the solid 
line is the asymptotic theory ([9|. 



Numerical simulations of the switching rate for purely 
Poisson noise are shown in Fig. [2l There is good agree- 
ment between the data and the asymptotic theory ([9]) for 
small g. In this range Rpg slowly varies with g. 



In conclusion, we have considered switching in systems 
simultaneously driven by a Gaussian and a non-Gaussian 
noise. Even where the non-Gaussian noise has intensity 
smaller than that of the Gaussian, it may strongly change 
the switching rate. The effect is determined by the ra- 
tio of the non-Gaussian noise amplitude to the Gaussian 
noise intensity. It is described by the characteristic func- 
tional of the non-Gaussian noise calculated for a function 
determined by the system dynamics without this noise. A 
non-Gaussian tail of the noise distribution may strongly 
modify the switching rate even for a small noise inten- 
sity. We demonstrate this effect using Poisson noise as 
an example. Analytical results are compared with Monte 
Carlo simulations. 
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